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ABSTRACT 


The general operation of KATE, an artificial intelligence controller, is 
outlined. A shuttle ECS demonstration system for KATE is explained. The 
knowledge base model for this system is derived. An experimental test pro- 
cedure is given to verify parameters in the model. 
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SUMMARY 


A scaled down version of the shuttle ECS system is being built to test KATE, 
an artificial intelligence expert system at KSC. KATE requires an accurate 
mathematical model of the ECS, called the knowledge base. This report gives 
the model derivation. 

An explanation of how KATE works is given. Each component must be described 
and loaded into KATE as a LISP "frame*. The description includes a func- 
tional equation and component relationships to aid fault diagnosis. 

The ECS system consists of four ducts branching from a manifold which is 
supplied with chilled air. Each duct has a heater, flow meter, a butterfly 
control valve, and two manual (butterfly) valves. The air temperature and 
flow rate is to be controlled precisely by KATE, Component failure is to be 
initiated to test KATE's diagnostics. 

The mathematical description relates flow rates to pressures. A major ele- 
ment is the determination of component loss coefficients. The heater loss 
coefficient and time constant were derived from the fundamentals of flow over 
fins. Butterfly valves show little control at angles less than 40° and too 
much control at angles greater than 60°. Manual valve settings in each duct 
are obtained which will optimize the valve control sensitivity. 

An experimental procedure is given to verify the loss coefficients of key 
components. Loss coefficients are expressed in terms of pressure and flow 
rate readings . 
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SYMBOLS AND SUBSCRIPTS 


Symbols 

a 

A 

C 

D 

f 

9r 

H 

K 

L 

M 

P 

T 

V 

a 

B 

P 

e 


admittance 
duct area 
constant 
duct diameter 

pipe friction factor, 0 f 017 

32.17 Ibm-f t/lbf - sec 2 

heat transfer rate 

loss coefficient 

duct length 

mass flow rate 

pressure 

temperature , Rankine 
velocity 

venturi loss, fraction of AP 
D 2 /D x 

air density 

butterfly valve setting or heater divergent wall angle to 
horizontal 



Subscripts 

a 

b 

c 

d 

e 

f 

h 

L 

m 

P 

r 

s 

T 

1 

2 


aft 

bend 

cabin 

dump 

entrance 

forward or fins 

manual (hand) valve 

line (duct) 

motorized (control) valve 
payl oad 

reheat (heater) 

supply 

total 

venturi entrance 
venturi throat 


w 
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I . INTRODUCTION 


KATE (Knowledge -based Autonomous Test Engineer) is a knowledge -based expert 
system being developed at Kennedy Space Center. The basic concept of an 
expert system consists of two parts. The first part is a knowledge base 
which describes the process hardware. This allows for a simulation of real 
time behavior , and it can be used to gather information about the status of 
the system. The second part is the software or artificial intelligence (AI). 
AI represents a good design engineer, systems engineer, and operator in a 
single software package that operates on the knowledge base representation of 
the system. Software capabilities include graphical display generation, sim- 
ulation, process control, redundancy management, constraint checking and 
diagnostics. In particular, unanticipated failures can be diagnosed, and 
instructions can be given to get the system to a new desired state, or cor- 
rect it and return it to a normal state. 

KATE has been under development at Kennedy Space Center since 1985. The 
first full test of the system will be carried out this Fall. The process 
chosen for the task is a scaled down version of the shuttle environmental 
control (ECS) ground system. Such a system consists of thermal, fluid and 
electrical components. In addition, the sensors and control actuators are 
both analog and digital. A mathematical description of these elements con- 
stitutes the process knowledge base. The summer task is to assist Kennedy 
engineers obtain this knowledge base. 



II. HOW THE SOFTWARE WORKS 


Artificial intelligence (AI) systems are composed of three parts: (1) 

input/output, (2) knowledge base, and (3) operating system. Input consists 
of the sensors (transducers), supplying the system with pressure, tempera- 
ture, flow, etc. data in the form of voltages to an analog -to -digital conver- 
ter. Outputs are the commands to actuators regulating the variables of 
temperature, flow, etc. The knowledge base is a mathematical description of 
the process. The operating system is the set of algorithms which interprets 
the data and knowledge base, controls the process and diagnoses faults. The 
latter is KATE. Her language is LISP although other AI languages include ADA 
and C. 

2 . 1 KNOWLEDGE BASE 

The operating system is sufficiently general that it can be applied to all 
processes. It is the knowledge base that changes with each process. There 
are two types: expert rules and model base. 

With expert rules, the programmer sits down with the plant operator and 
grills him on how he controls the process. A plant is usually modified over 
the years, and the operator may be the only person who knows the system. He 
is aware of the system idiosyncrasies and often learns by experience why a 
certain procedure works. The programmer must root out these rules of intu- 
itive behavior. These rules are of the form: if , then . For 

example, if valve one is open or valve two is closed and liquid level switch 
two is set, then open valve three. The rules are adaptable to Boolean alge- 
bra. Some large systems have been known to take five man-years and generate 
thousands of rules. Drawbacks to this approach are several. 

a. It cannot react to unknown conditions. 

b. There is no diagnosis when faults occur. The process must be shut 
down with a simple sensor failure. 

c. The intelligence and system operation resides mainly within the 
knowledge base rules; and therefore, they change with each process. 

KATE uses a "model base" knowledge base. It depends on a mathematical 
description of each component in the process. Thus, each measured variable 
can be calculated also. A comparison of the two is the basis for a complete 
fault diagnosis. At the same time the AI system can suggest alternative 
commands to circumvent the failed component. Intelligence resides in the 
operating system, and only the component mathematical descriptions need to be 
loaded when a new process is brought under AI control. 
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As an example of a component description, consider a venturi flow meter in an 
air duct. Normally, one calculates the volume flow rate from the differen- 
tial pressure reading. In our case we want the mass flow rate which will 
include the duct temperature and pressure (in place of density) calculation. 
The equation may look like 


DP -6.2PM 7 /T 

Notice that the equation is the inverse of the normal output M, i.e. we want 
the differential pressure in terms of the flow rate. The reason is that the 
AI system must calculate DP to compare it with the measured DP for fault 
diagnosis . 

All components and variables are prepared for programming by placing its 
information into a "frame". An example follows for the differential pressure 
above . 


LISP FRAME 


(deframe DPA 

(nomenclature aft duct differential pressure) 
(a i o PRESSURE) 

(source -path P 
T 
F) 

(in-path-of DPA TRANSDUCER) 

(status (/(*(6.2)( square Ai)(P))T))) 


Line 1. Variable symbol. 

Line 2. Description. 

Line 3. An instance of. What is the category? A variable: pressure, 

temperature, flow? An item that can fail: transducer, relay, 

heater? 

Line 4. Source -path. What variables are used to calculate DPA? 

Line 5. In-path-of. What component is in the path of DPA? 

Line 6. Status. The equation for DPA. 

2.2 KATE'S OPERATING SYSTEM 

KATE'S operating system consists of algorithms for fault diagnosis which 
function basically by forming lists. After a delay for command dynamics to 
settle out, the following procedure is used to identify failed components and 
determine alternate commands to circumvent the problem. 


a. Get new measurements from sensors throughout the process. Compare 

these new measurements with previous readings stored in memory. Form 
a list of any measurements that change by more than a prescribed 
amount . 
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b. Go to the knowledge base and calculate the variable that is measured 
for each discrepancy in the above list. If the calculated and mea- 
sured values disagree, place the deviants into a second list for the 
DIAGNOSER. 

c. When a component fails, it is likely to generate many deviants for 
the second list above. The DIAGNOSER goes to the knowledge base and 
forms a list of all possible components that can be related to the 
deviant. Recall the frame for differential pressure. Line 5 identi- 
fies a transducer as a possibility. Line 4 has several variables 
that can be related to components upstream through other frames. The 
computer knows which categories (Line 3) in every frame that can 
fail . Every deviant measurement caused by the failed component will 
have a similar, though not identical, list of process components that 
factor into its reading. Every list will have several process compo- 
nents in common., KATE uses elementary set theory to identify those 
components with the most intersections and ranks them in the order of 
most likely to cause the deviations. Starting with the first compo- 
nent in the list, KATE simulates the component failure modes and 
calculates every measurement in the process from the knowledge base. 
The component mode can be ON/OFF for a relay to multiple command 
voltages for a motor control valve. KATE proceeds through the list 
until a component and its failure mode is found which will cause 
agreement between all process measured and calculated variables. If 
KATE is told to find a new command to circumvent the failure and to 
maintain variable settings, the procedure is the same as above. Com- 
mand changes are simulated until the desired result is obtained * 

This is valuable at Kennedy Space Center where systems often have 
built in component redundancy. 
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III. SHUTTLE ECS DEMONSTRATION UNIT 

A top view drawing of the scaled down shuttle environmental control system 
(ECS) is shown in Figure 1. It is built within a 18 by 13 foot open block- 
house. An external chiller (a mobile air force purge unit) supplies 
100 lbm/min of cold air at 43° F, 3 to 4 psig through an 8 inch line. The 
line turns 180° into a 12 inch manifold which distributes the cold air to four 
ducts for the simulated shuttle payload, aft, forward and cabin areas. A 
motorized dump valve balances the out flow for the total of 100 lbm/min. 

Each duct is identical and has the same components with one exception: the 

payload duct is a nominal 6 inch pipe (6.065ID) while the remaining ducts are 
a nominal 3 inch pipe (3.068ID). The components are as follows. 

a. Keystone hand operated butterfly valve. 

b. Chromalox air heater (6kw and 15kw, payload duct). 

c. Leeds and Northrup venturi meter. 

d. Keystone motorized butterfly valve. 

e. Keystone hand operated butterfly valve. 

Each duct also contains, as indicated in Figure 1, the following. 


a. Two remote pressure sensors. 

b. One remote temperature sensor. 

c. Two pressure gauges. 

d . One temperature gauge . 

e. One remote differential pressure sensor (venturi). 

Both the air temperature and flow rate are controlled in each of the four 
ducts. The temperature is controlled to within + 0.5° F over the range, 
to 70° F. The flow rate is controlled to within 0.2 lbm/min where the range 
for each duct is as follows. 


Payload duct: 
Cabin duct: 
Forward duct : 
Aft duct: 


28 to 55 lbm/min 
4.6 to 12 lbm/min 
6 to 17.6 lbm/min 
11.4 to 14.2 lbm/min 
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IV. ANALYTICAL DESCRIPTION 


4.1 BASIC EQUATIONS 


A mathematical description of the overall test unit consists of equations 
relating the flow rate and pressure drops for each duct in terms of admit- 
tance . 


M.-aJP.-P., 


M p “ a p'/P* 
M,-a,/P^ 


M'-ct'JK 


M s -a,JP~ n 

a,-a p + a a + a / + a c + a < | 


Each duct includes a number of components and elements in series 
admittance can be found in terms of the individual admittances . 
the aft duct admittance is found from 


The duct 
For example, 


a 


a 



1/ 



)/? 


vhere n is an individual admittance . 

The admittances for each element can be expressed in terms of a loss coeffi 
cient K. The pressure drop is given by 

AP • KpV z /2g c 
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where M - pAV. Rearranging, 


a = A^2pg^K 


For a perfect gas, 


P/Po 


LTi 

r, t 


(use absolute values) 


Then 


a - y[2p , 9 c /I fFTjTj / 


For air at 70° F, 14.7 psia, 

a - 1 1 .0/1 v'PT./P.T / /K lbm/min /psi 
a = 2.09/1 f?T c /P Q T / lbm/min ,/in .« 2 0 
and for elements in series 


a - ll.Od/PTvTj 1 [£*,]' 1/2 

We will now review the individual components and elements in order, starting 
with the purge unit. The aft duct branch will be taken but the calculated 
variables are the same for each duct unless stated otherwise. 

4.1.1 SUPPLY PIPE. 


K sl ~ }L/D- (0.01 7)(20) / ( .666) - 0.5 
4.1.2 SUPPLY FLOW METER. 


/ 2pg c A/ > 7 

M,-YC„yi 2 J - B «- (AP t across taps) 

Y (compressibility factor) * 0.98 (neglect) 
AP„ (loss across venturi) - aAP T 

a = 0.10' for 0 = 7°, B = 0.45 
Flow Measurement Engineering Handbook REF TC 177 
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C d -0.96 



Loss between P, and RP9, 


Y.K-Ksi + K s* m 3 ° 


4.1.3 90® BEND. 


*s." 0.3 


4.1.4 EXPANSION TO MANIFOLD. 


K 



.3 


4.1.5 CONTRACTION: MANIFOLD TO DUCT. 

“ 0.45 

-0.36 (PLB) 

4.1.6 BUTTERFLY VALVE (MANUAL 1). 

[ 1-56 V 

Ao/ii “ l j _ s inO ) 


e>25® 


See Section 3.1.12 

4.1.7 REHEAT CHAMBER. 
Gradual enlargement 



.9 


Fin drag 


1.328 



x 400 fins 
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Cylinder drag 


pVl 

F ' mC * 2^ c A <* 6 Cyh 

Loss coefficient due to drag 

A ?A k -(F,+ F c ) 

also 

&P- K 2 pVl/2g 


K 2 


A / 


1^400. c a A c x 6 

4*1 


- 76 


correct to duct velocity at A 

K z - 76( A r / A) z m .25 

total reheat chamber loss 

K ar - o.9 + 0.25= 1.2 (not calculated for PLB but scales up similarly) 

reheat chamber time response 

H MCdT ,/di heat transfer to fin capacitance 

H 2 ~(T f -T 0 )h / A / heat transfer to film resistance 
H - H , + H 2 - MCdl ,/ dt + h f A t {T f - I B ) 

But heat transfer to the film resistance is given to the air. 

( T r T^h f A f -M a C p {T t -T , ) 

Eliminating T , and rearranging H, 

A f C tt ( 1 + (M a C p / h f A d (T 0 - T ,) / dt + M a C p (T 0 - T J - fl 
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The time delay t d is four time constants or 


*d 


4 MC 
M a C p 


(l*(M a C p /h,A,))-3.25 


min 


-4.0 (50% fin efficiency) 


where 


h, - 0.664 j/ > i /3 (*.i) 1 ' 2 - 4* 10 

^y-840 in z (est.) 
C-0.11 Btu/lbm°F 
M-ll.b lbm (est.) 
C p - 0.24 


reheat chamber sensitivity 


Power - M a C,(T 0 - 7,) 

at full power the temperature change of the air is 

60001/ x 3.414 ^ , n . 

T °~ T, ‘\? lbm/min x 0.?4 Btii/lbrnT 

For an 8 bit D/A converter, the temperature can be controlled to within 0.5 F 
at flow rate mid-range. Worse case: low cabin flow rate (4.6 lbm/min) where 

T 0 -t,“ 308" or control to within 1.2° F. 

4.1.8 DUCT SECTION HEATER TO RP4 SENSOR. 

K ell -fL/D-0.18 

-0.09 (PLB) 


Loss between RP9 and RP4 , 


£ K - K Sl + K „ +* + K ohl * K or * 


- 2.4 + K ah) 

« 

* 
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4.1.9 DUCT SECTION SENSOR TO FLOW METER. 


K al ,-fL/d- 0.1 
-0.05 (PLB) 


4.1.10 FLOW METER. 

Sane equations as the supply flow meter. 

3 in. duct : B-0.671, 0-5" 

6 in. duct : B-0.4G, 0-7° 


Meters for ducts 4 inches in diameter or less have a slight sudden contrac- 
tion from duct size A to area A,, at the high pressure tap. 


K 




1.6 


K 2 -.28 (sudden contraction) 


K av -K^K Z - 1.9 

^-^( p )’ 25 (PLB) 

4.1.11 DUCT SECTION TO MOTORIZED BUTTERFLY VALVE. 

K al3 -fL/D- 0.1 

-0.05 (PLB) 


4.1.12 MOTORIZED BUTTERFLY VALVE 

If both hand butterfly valves are full open, a simple calculation will show 
that the motorized valve must be set to nearly closed 0*70" for the mid-flow 
range where its A." * 600. A one degree change will change K by 150, much too 
sensitive for satisfactory flow rate control. Furthermore, at nearly open 
0-20°, a one degree change will change K by 0.5. If the hand valves take too 
much of the loss to achieve the mid-flow rate setting, the motorized valve 
loses control . A logical solution is to select a minimum permissible setting 
for the maximum flow rate. This will maintain control, yet allow a maximum 
sensitivity for control. Selecting G min “40° 


11 
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No. we can determine the total * by summing the Ks along a streamline fro. P. 
Spurge unit, to a selected duct end at r-Opslg. The Ks In the supply 
duct should be corrected to the common velocity of a distributing 


K.,d* 






Then 


M a -aJF, 


.here o - 11 .8*1* using an average value for the pressur./temperature 

correction. The butterfly setting for the two hand valves Is found fro. 
maximum duct flow with the control valve at 0 - 40 . 


f MV 


11.84120 + 2 *,]^. 


solving for each duct (/>,- 3.5 psig) 


50 

0-53.8' 

K^-K^-29 

0-49.1' 

*t»l “ Kck 2 " 75 

0-57’ 

* phi m % pA2 “ 

0-54.4* 


Butterfly control valve settings for the aver.ge flo. rat. In each duct (hand 
valve settings above) are 


K a ,-40 

0-51.9’ 

K /m -108 

0-59.7' 

^ci-1 94 

0-63.6' 

* pm “ 1 09 

0-59.7’ 

4.1.13 DUCT SECTION CONTROL VALVE TO 

SENSOR RP3. 

^ ota “ 

/I/D- 0.10 
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4.1.14 DUCT SECTION SENSOR RP3 TO HAND VALVE. 

* oi4 - /I/D- 0.10 

4.1.15 BUTTERFLY VALVE (MANUAL 2). 

f 1.56 y 
K °™~{ l -sinG ) 


See Section 3.1.12 


0 > 25" 
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V . EXPERIMENTAL DESCRIPTION 

An accurate model of the process begins with an analysis of the components 
and a calculation of the most probable values for the variables. Along the 
way certain engineering parameters must be estimated. The final link in the 
model development is an experimental program to verify the analysis and its 
associated parameters. For the ECS system there are several verifications 
that are most important. 

(1) Characteristic curve of the purge unit: /‘.-/(M,) 

(2) Loss coefficient versus angle for the butterfly valves. 

(3) Flow meter loss coefficient. 

(4) Heater loss coefficient. 

(5) Heater response time. 

The ECS system is now being assembled and it should be completed for tests by 
the middle of August. The initial tests to verify the parameters will be run 
without KATE, using the various gauges that accompany the remote sensor 
units. I propose the following test procedure while recording all gauges. 

5.1 OPEN ALL DUCT BUTTERFLY VALVES 

a. Record minimum P s on the characteristic curve along with its M,. 

b. Find the open manual butter fly loss coefficient K ahZ . It should be 
the same for all valves. 


M a -afPG3 


a-c[A' oM + A 0 , 2 ( open)] 


c. Find the venturi loss coefficient K ov 


M 0 “ a PG 4 _ PC 3 

a - c[K aL2 * K av * K oL3 * *„„( open) + K ol3 ] 


289 



d. Find the heater loss coefficient K ar 

m a - a^/ > G9- PG4 

<fc[(K SB + K„)] .4 + K ahl {ope.n)+ K ar + K all Y uz 

e. Check the supply admittance. 

Ai, - a^/P,- PG9 

a-c'[K sn + K„+K. I2 r , ' z 

f. Check all ducts for identical gauge readings. 

5.2 VARY CONTROL VALVE SETTINGS IN CABIN DUCT. 

a. Obtain A c „-/(0) 

b. At A cm (closed) check to see if manifold pressure given by PG8 agrees 
with P s . This verifies that the supply line loss Is negligible. 
Record P s , M s for another point on the characteristic curve. 

5.3 CLOSE REMAINING DUCTS IN TURN 

a. Completely close aft, forward and payload bay ducts, recording P s and 
M s after each closing. 

b. Slowly close the dump valve. Continue recording P, , for the 
characteristic curve. Stop when P , reaches a maximum. 

5.4 DETERMINE HEATER TIME RESPONSE 

Introduce a step change in the heater command. Measure the outlet air tem- 
perature every 10 seconds until a new temperature stabilizes. 



w 
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VI. CONCLUDING REMARKS 


KATE depends on an accurate knowledge base flow model . The ECS system is 

SpL with various loss coefficients the predo.in.nt determinate 
Among the loss coefficients calculations for hends. cont , 1 
s.ooth pipe lengths and even venturi Jlo. .et.rs are str. * these 

handbooks. These values are small and very rename. y f or 

contributions between pressure gauges, the less certain loss co heater 

the heater and butterfly valves can be verified experl.ent.lly The heater 
loss coefficient was calculated fro. aerodyna»lc drag principles. The fo 

IS! for the butterfly valve loss coefficient versus angle was found^in^he 

■Handbook of Hydraulic Resistance- . translated fro. R “»» l • 

indication of the formula source, likely to be theoretical. 

The single critical element in the knowledge base is the butterfly val ^«- 
Se valve angle is set by activating the butterfly drive motor for a speci- 
fied time. The multiple spur/worm gear drive is iMy position 

Thus open loop positioning cannot be reliable, kaie. snouxu r 

potentiometer Suh some type of software feedback . 

bv the butterfly flow sensitivity. The mid-range flow rate valve settings 
for the three series (two manual plus one control) i butterfly va are 
approximately 55°. These settings were established to give e 
allowable change in ioss coefficient^ angle rf.ove.ent, 

^high'fS Si will change°by 0,! 5 lb./.in per degree while ^ low flow 
rates will change by 0.55 lbm/mln per degree over the angle span of 30* (for 
ward duct example). 

Finally the heater time response was calculated to be epproxiMtely 4 min- 
utes independent of the temperature command change magnitude. The length o 
“S Sbe large. ?hus, it represents an element of anticipation in 

the experimental study . 
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